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ABSTRACT 

An  adaptive  control  system  operating  on  submarine  depth  rate 
with  constraints  on  the  maximum  values  of  stern  plane  deflection 
angle  and  hull  pitch  angle  is  formulated  using  a  linear  compensation 
technique.   The  system  is  adaptive  in  the  sense  that  it  is  designed 
to  give  near  uniform  dynamic  response  over  ranges  of  values  of  various 
hull  geometry  and  speed  parameters.   The  formulated  system  is  simulated 
on  a  Reeves  Electronic  Analog  Computer,  various  adjustments  made  and 
justified,  and  results  shown  in  plots  of  stern  plane  angle,  pitch  angle 
and  depth  rate  error.   A  simple  system  for  predicting  pull-out  time  is 
suggested  which,  when  associated  with  the  proposed  system,  would 
facilitate  rapid,  well-damped  depth  changes. 
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Symbols  and  Nomenclature 


The  nomenclature  defined  in  T  &  R  Bulletin  No.  1-5  of  the  Society 
of  Naval  Architects  and  Marine  Engineers  is  used  herein  where  applic- 
able.  A  right-handed  cartesian  coordinate  system  is  employed  with  the 
origin  at  the  center  of  gravity  of  the  hull,  the  positive  x  axis  dir- 
ected forward,  the  positive  Z  axis  directed  down  normal  to  the  base 
line,  and  angles  measured  positive  in  the  counter-clockwise  direction. 
Pitch  angle  and  depth  are  measured  relative  to  a  similar  system  fixed 
in  space  and  oriented  with  the  positive  Z  axis  directed  vertically 
down.   Quantities  measured  relative  to  this  fixed  system  are  designated 
by  the  subscript  g.   All  quantities  are  non-dimensionalized,  as  indi- 
cated by  the  prime,  by  dividing  by  appropriate  powers  of  the  length  of 
the  hull  and  linear  velocity,  the  mass  and  mass  density. 
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Z-j  Z.J,    -   metacentric  height 


r  S    '  &        rad  stern  plane  angle 
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i  desired  or  commanded  quantity 

o  actual  or  output  quantity 
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INTRODUCTION 

A.   General  Aspects  of  the  Depth  Control  Problem 

Prior  to  the  closing  years  of  World  War  II,  submarine  design  features 
were  based  primarily  on  surface  performance.   As  a  direct  result  of  this 
philosophy,  submerged  speeds,  and  hence  maneuverability,  were  slow  and 
sluggish.   The  top  speed  of  a  typical  boat  of  this  period  was  only  nine 
knots  at  the  half-hour  battery  rate.   For  this  reason  the  crude  two- 
plane  sman  depth  control  system  generally  employed  since  the  very  earliest 
days  of  submarine  construction  was  entirely  adequate. 

The  improved  methods  of  anti-submarine  warfare  brought  about  by 
developments  in  sonar  and  radar,  the  use  of  carrier-bourne  aircraft, 
and  the  construction  of  fast  anti-submarine  escort  ships  made  higher 
submerged  speeds  imperative.   Belatedly,  submerged  performance  became 
the  primary  consideration.   By  greatly  reducing  the  amount  of  super- 
structure, and  fairing  the  bridge  and  periscope  sheers,  hull  resistance 
submerged  was  reduced  by  almost  a  factor  of  two.   This,  coupled  with  a 
doubling  of  battery  capacity,  increased  submerged  speed  by  a  factor  of 
two.   These  higher  speeds  indioated  clearly  the  need  for  faster,  more 
accurate  control  of  depth.   This  need  is  easily  recognized  when  one 
considers  that  submarines  of  this  period  would  collapse  at  depths  little 
in  excess  of  one  ship  length.   Even  present-day  submarines  with  their 
much  higher  speeds  must  confine  their  motions  to  a  layer  of  water  only 
a  few  ship-lengths  deep. 

The  advent  of  nuclear  propulsion  plants  and  the  low  resistance, 


body  of  revolution  hull  form  developed  in  the  last  decade  have  greatly- 
magnified  this  depth  control  problem.   Let  us  consider  the  situation  con- 
fronting the  submarine  pilot  controlling  the  stern  planes  when  only  stern 
planes  are  used  to  control  pitch  angle  and  depth.   This  situation  is 
adequately  represented  by  the  following  blook  diagram. 


Zgi 

commanded  depth 

Ss 

stern  plane  angle 

Zgo 

actual  depth 

Pig.  I 

e 

pitch  angle 

From  this  diagram  it  is  evident  that  there  are  three  integrations  be- 
tween the  immediate  effect  of  the  operator's  action  and  the  desired 
result.   If  the  operator  is  to  close  the  control  loop  in  an  optimum 
manner,  he  must  anticipate  the  results  of  these  integrations  -  a  matter 
difficult  in  itself,  but  greatly  complicated  by  the  long  lag  between  0 
and  §  introduced  by  the  la'ge  hydrodynamic  mass  and  small  hydrodynamic 
stiffness  of  the  hull.   The  operator  attempts  to  overcome  this  diffi- 
culty by  introducing  a  large  initial  §    .   This  gives  him  a  more  rapid 

s 

indication  of  the  result  of  his  control  action,  but  due  to  the  lag 
mentioned  above,  he  has  no  timely  indication  of  when  to  remove  or  re- 
verse plane  angle.   This  results  in  overshooting  the  desired  objective 
-  the  behavior  of  a  poorly  damped  system.   Thus  it  is  evident  that  the 
operator  is  required  to  perform  two  functions  simultaneously  in  closing 
the  control  loop:  he  must  speed  up  the  system  response,  and  he  must 


provide  adequate  damping  characteristics.   This  latter  skill  requires 
a  certain  amount  of  prescience  highly  dependent  not  only  on  the  skill 
of  the  operator,  but  also  on  his  physical  condition.   In  this  sense, 
operator  reaction  is  a  random  variable  introduced  into  the  control  pro- 
blem which  precludes  absolutely  consistent  performance. 

Two  basic  approaches  have  been  followed  to  alleviate  or  avoid  these 
difficulties  in  depth  control,,   The  first,  and  most  apparent,  is  replace- 
ment of  the  human  operator  in  the  control  loop  by  an  automatic  plane 
angle  computer.   The  second  approaoh,  and  the  one  followed  in  this  in- 
vestigation, involves  relocation  of  the  human  operator  outside  of  the 
closed-loop  portion  of  the  system. 
B.   Development  of  Automatic  Depth  Control  Systems 

The  German  Navy  was  the  first  to  practice  the  new  philosophy  of 
submarine  design.   Their  type  XXI  boats  constructed  just  before  the  end 
of  World  War  II  in  Europe  were  capable  of  about  twenty  knots  submerged. 
These  ships  clearly  indicated  the  desirability  of  an  automatic  depth 
control  system.   An  attempt  was  made  to  meet  this  need.   The  system 
which  evolved  was  really  two  very  similar  sub-systems  -  one  for  each 
set  of  planes.   Plane  angles  were  computed  as  functions  of  speed, 
depth,  and  pitch  angle.   Depth  and  pitch  angle  signals  were  derived 
from  pressure  gauges  and  pendulums  respectively.   Variable  gain  settings 
in  each  of  the  sub-systems  permitted  ratioing  bow  and  stern  plane  angles 
as  a  function  of  speed  and  depth-keeping  accuracy  requirements.   The 
controlled  quantity  was  depth,  so  that  the  overall  system  was  basically 
a  depth  regulator  rather  than  an  autopilot  capable  of  radical  maneuver- 
ing. 


The  United  States,  profiting  from  German  experience,  began  modifying 
the  top-hamper  and  battery  capacity  of  their  SS-475  class  fleet  boats  in 
the  late  1940* s.   Coincident  with  this  Guppy  conversion,  development  of 
a  semi-automatic  depth-keeping  system  was  initiated.   The  first  such 
system  was  built  by  Askania  Regulator  Company  and  installed  on  the 
U.S.S.  Irex  (SS-482)  in  1948.   Although  very  similar  to  the  German 
system  in  most  respects,  it  controlled  not  only  depth,  but  also  boat 

angle  or  pitch  angle.   With  this  added  feature,  it  provided  more  man- 

5 
euvering  capability  than  the  German  system.   The  first  of  these  systems 

were  semi-automatic  in  that  the  output  was  an  indicated  plane  angle 
which  the  planesman  was  required  to  match.   Subsequent  systems  were  pro- 
vided with  an  automatic  follow-up  device  making  them  fully  automatic. 

With  the  advent  of  Albacore  and  submerged  speeds  in  excess  of 
twenty  knots,  these  early  systems  became  inadequate.   Present  systems 
are  fully  automatic,  and  similar  to  aircraft  autopilots.   The  Sperry 
System,  installed  on  most  new  construction,  determines  stern  and  sail 
plane  angles  from  signals  in  depth,  depth  rate,  pitch  angle,  pitch  rate, 
and  pitch  acceleration  to  keep  the  boat  on  a  digitally  computed  optimum 
trajectory. 
C.    Developments  in  Manual  Control  Systems 

The  high  submerged  speeds  and  limited  depth  range  of  modern  sub- 
marines make  the  poor  coordination  and  large  time  lags  inherent  in  the 

*The  Guppy  (Greater  Underwater  Propulsive  Power)  conversion  involved 
streamlining  the  superstructure  and  doubling  the  battery  capacity 
of  the  SS-475  class  boats.   These  two  modifications  resulted  in  a 
200$  increase  in  submerged  speed.   The  resulting  ships  bear  an  un- 
canny resemblance  to  the  type  XXI  boats. 


old  two-operator,  manual  control  system  intolerable.   The  first  step 
away  from  the  old  system  was  incorporated  in  Albacore  with  the  substi- 
tution of  an  aircraft- type  control  column  for  the  old  wheel  type  de- 
vices.    Two  sticks  were  providedo   The  main  unit  is  a  conventional 
aircraft  steering  column  in  which  wheel  rotation  controls  rudder  posi- 
tion, and  fore  and  aft  motion  of  the  entire  column  controls  bow  plane 
angle 0   The  second  stick  controls  stern  plane  angle.   In  the  principal 
mode  of  operation,  stern  plane  angle  command  is  derived  from  the  bow 
plane  signal  via  a  variable  gain  amplifier,  thus  providing  single- 
operator  ratioed  control.   This  is  the  mode  of  operation  greatly  pre- 
ferred by  the  operatorsc 

This  basic  system  has  proven  satisfactory  in  most  respects.   Sub- 
sequent developments  have  been  primarily  concerned  with  providing  the 
operator  with  the  required  feed-back  information  -  depth  error  and  boat 
angle  -  in  the  most  efficient  manner.   The  first  major  improvement  in 
this  line  was  the  Combined  Instrument  Panel  which  displays  all  the 
information  required  by  the  operator  in  a  cluster  of  dials  arranged 

Q   -I  -| 

in  an  optimum  manner  from  the  human  engineering  point  of  view.  '    A 

second,  and  rather  novel,  development  is  the  Contact  Analog  Display 

o 

developed  by  the  Electrio  Boat  Division  of  General  Dynamics,,   In  es- 
sence, this  is  a  television  screen  picture  of  the  path  the  operator  is 
required  to  follow.   A  three  dimensional  effect  is  provided  by  a  fic- 
ticious horizon  and  narrowing-with-di stance  perspective  of  the  pic- 
tured road.   This  is  the  only  indicator  the  pilot  must  concentrate  his 
attention  upon.   His  task  is  simply  to  keep  the  white  line  of  the  road 
properly  positioned  on  the  screen. 

In  anticipation  of  the  performance  capabilities  to  be  required  of 


future  high  speed  submarines,  the  Bureau  of  Ships  has  initiated  an  ex- 
tensive study  of  the  over-all  submarine  control  problem.   The  title  of 
this  study  is  the  Submarine  Integrated  Control  Program  (SUBIC).   The 
Contact  Analog  mentioned  above  originated  in  this  program.   As  a  part 
of  this  investigation,  Cornell  Aeronautical  Laboratories  conducted  a 

study  of  the  data  flow  and  data  processing  requirements  in  effective 

7 
submarine  control.   One  of  the  major  results  of  this  study  was  an 

indication  of  the  desirability  of  locating  the  human  operator  external 
to  a  multiple  closed-loop  system  in  which  the  outer  loop  closure  is  on 
depth  rate.   The  advantages  of  this  scheme  are  readily  apparent.   In 
addition  to  minimizing  the  effects  of  erratic  human  behavior  and  sim- 
plifying the  operators'  problem  by  avoiding  two  of  the  integrations 
mentioned  in  (a)  above,  the  system  is  easily  made  fully  automatic  by 
substituting  a  digital  computer  or  programmer  in  place  of  the  human 
operator.   This  system  is  represented  by  the  following  block  diagram. 
The  K's  indicated  are  in  general  non-linear  functions  of  ship  speed  and 


— * 

Pilot 

p  . 

1L 

— * 

Computer 

KXJ 


Fig.  II 
the  geometrical  parameters  of  the  hull. 


I).    The  System  to  be  Investigated 

The  major  difficulty  in  submarine  control  system  design  is  the  lack 
of  precise  information  concerning  the  values  of  the  various  hydroframe 
parameters.   At  the  present  time,  there  is  no  satisfactory  method  of  es- 
timating these  parameter  values  within  acceptable  accuracy  from  the 

geometry  of  the  hull  as  delineated  on  the  lines  drawing.   Recourse  to 

12 
model  and  full-scale  tests  is  unavoidable.    Further,  in  order  to 

obtain  acceptable  performance  over  the  speed  range  of  a  given  hull,  the 

characteristics  of  a  control  system  such  as  that  shown  in  Fig.  II  must 

be  varied  in  some  pre-determined,  generally  non-linear  manner  with 

speed.   It  is  with  this  problem  that  this  thesis  is  concerned. 

The  basic  system  under  consideration  is  shown  in  Fig.  III.   It  is 

evident  from  this  diagram  that  no  provision  is  made  for  a  steady  state 


Fig.  Ill 
plane  angle  in  a  steady  dive.   This  feature  would  require  an  auxiliary 
forward  loop  from  the  input  to  the  plane  angle  summer.   This  matter  is 
ignored  in  the  following  investigation  because  the  effect  of  such  a  for- 
ward loop  is  small  and  decreases  with  speed  roughly  as  the  reciprocal  of 
the  speed  cubed.   The  effects  of  this  simplification  will  be  a  slight 
increase  in  steady  state  depth  rate  error  at  slow  speeds.   Also,  only 


stern  plane  control  is  considered. 

The  objective  of  this  investigation  is  to  formulate  suitable  com- 
pensation functions  K  .G  n  and  K  ^G  _  that  will  yield  a  specified  per- 

cl  cl      c2  c2  *        *  * 

formance  of  the  closed  loop  system  over  a  specified  speed  range  of  a 
given  hull  configuration,  and  at  the  same  time  minimize  the  effects  of 
deviations  of  certain  important  hull  parameters  from  their  specified 
values.   In  this  sense,  the  proposed  system  is  an  adaptive  system. 

The  problem  will  be  treated  on  a  non-dimensional  basis  using  the 
linearized  force  and  moment  equations  of  a  submerged  body  in  the  ver- 
tical plane,  and  standard  nomenclature  as  specified  in  reference  \_1~]. 
The  system  will  be  simulated  on  a  REAC  analog  computer  to  determine  its 
transient  and  steady  state  response  to  a  step  input  of  depth  rate  under 
various  conditions. 


II 

PROCEDURE 

A.    Scope 

This  investigation  is  concerned  with  the  formulation  of  a  suitable 
closed-loop  system  for  controlling  the  depth  rate  of  a  deeply  submerged 
submarine  subject  to  constraints  on  the  maximum  values  of  hull  pitch 
angle  and  plane  deflection  angle.   The  principal  objective  is  a  system 
whose  dynamic  response  to  an  ordered  change  of  depth  rate  is  relatively 
independent  of  variations  in  ship  speed  and  hull  geometry.   A  linear 
design  technique  will  be  employed,  and  the  performance  of  the  resulting 
system  analyzed  on  an  electronic  analog  computer. 

The  basic  idea  underlying  the  approach  to  be  followed  is  that 
hydroframe  parameters  may  be  considered  to  fluctuate  from  hull  to  hull 
as  well  as  for  a  given  hull  about  certain  normal  values.   Similarly, 
ship  speed  may  be  considered  as  a  deviation  from  some  normal  speed. 
The  effects  of  these  deviations  is  a  change  in  the  resonant  frequency 
of  the  hull  and  hence  of  the  entire  control  system.   If  when  its  para- 
meters have  their  so-called  normal  values,  the  overall  system  is  so 
designed  that  its  resonant  frequency  is  much  higher  than  any  frequency 
to  which  it  is  normally  required  to  respond,  then  the  change  in  natural 
frequency  brought  about  by  parameter  variations  will  have  a  small  effect 
on  system  performance. 

The  design  specifications  for  the  proposed  system  are  as  follows: 
1.  a  non-dimensional  natural  frequency  of  the  "normal"  system 
in  the  vicinity  of  5  with  a  damping  ratio  of  about  .75. 


2.  no  more  than  15$  overshoot  in  pitch  angle.* 

3.  no  more  than  10$  overshoot  in  depth  rate. 

4.  a  rise  time  to  10  pitch  angle  for  large  changes  of  depth 
rate  close  to  1  second  in  non-dimensional  time. 

The  non-dimensional  form  of  the  equations  of  motion  of  the  hull 
are  employed  since  this  facilitates  separation  of  the  dominant  para- 
meters from  those  of  lesser  importance. 
B.   Equations  of  Motion  of  the  Hull 


Omitting  the  cross-acceleration  terms  M.  w'  and  Z. ■ . ,  which  are 
°  w         q  q' 

2 
usually  quite  small,   the  linearized  equations  of  motion  of  a  submerged 

body  in  the  vertical  plane  may  be  written: 

(1)  (m'-Z.')  w'-  Z'  w1  -  (m'  +  Z' )  q'  -  Zl  S 
v/v    w  '     w     x     q      §s 

S 

(2)  -M'  w'  +  (I'  -  M'  )  4'  -  M'  q'  -  s£    ZI  6  =  M'  § 

v  '    w      v  y    q'       q  u    **   B      Sa      s 

s 

The  coefficients  of  the  acceleration  terms  w'  and  q'  are  inter- 
related through  the  following  equations: 


2 


m 


I'  -  (k«  ) 

y   y 

Z.'»  -k_m' 
w     2 

M, '  m   -k'  I' 

4       y 

k'  is  simply  the  non-dimensional  radius  of  gyration  of  the  hull  about 

the  y  axis.  (Cp  and  k'  represent  the  added  mass  of  the  hull.   They  are 

usually  taken  as  the  theoretical  values  applying  to  the  ellipsoid  of 

revolution  approximated  by  the  hull  form. 

*It  should  be  noted  that  the  desired  glide  angle  is  implicitly  included 
in  the  commanded  depth  rate  since  2  =  w  -  u60  Assuming  a  steady  glide 
with  zero  normal  force,  w  ■  o  and  Z°  =  -  U  6.  In  non-dimensional  nota- 
tion, Z'  =  -  6  8 

9    s 
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Landweber  and  Johnson  of  the  David  Taylor  Model  Basin  have  pro- 
posed the  following  relations  among  the  coefficients  of  the  velocity 

terms  for  the  case  under  consideration  where  the  only  appendage  is  a 

12 
tail  surface. 

M'   -  m'  -  X'  Z« 

w  1  w 

zf  =  x;  z« 

q      2   w 

M"   -   (X!)2  Z' 
q       3'   w 

These  relations,  although  only  approximate,  give  predictions  of  the 
indicated  parameters  which  are  reasonably  close  to  measured  values. 
x1  may  be  interpreted  as  the  non-dimensional  effective  moment  arm  at 
which  the  lift  due  to  angle  of  attack  acts.   The  m'  in  this  relation 
is  taken  as  an  approximation  to  the  Munk  moment  which  appears  in  the 
potential  flow  about  the  body.   X'   may  be  regarded  as  the  non-dimen- 
sional effective  radius  through  which  the  angular  velocity  q  is  related 
to  an  equivalent  normal  velocity,  i.e.  w1  =  X'  q' .   (X')   may  be  inter- 
preted as  the  effective  product  of  the  radius  relating  angular  and 
linear  velocities  and  the  moment  arm  relating  moment  and  force  due  to 
angle  of  attack.   Although  such  a  product  is  not  necessarily  a  square, 
it  seems  nearly  so  and  is  so  defined  for  convenience. 

The  coefficient  of  6  in  equation  (2)  may  be  reduced  to  a  speed 
parameter  as  follows:   m<?  Z,, 


-  a'  V 

—            — 

1            g      g 

$f>L  U     L 

an 

-      «ZB 

■ 

1 

m' 

rrr) 

U2 

~2 
U 

+M* 

This  speed  parameter  U  is  then  dependent  not  only  on  the  ship  speed  U, 
but  also  on  the  metacentric  height  Z_,  which  is  a  measure  of  the  ship's 
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longitudinal  static  stability.   Incorporation  of  the  minus  sign  in  the 
definition  of  l/u  is  done  since  Z_  is  always  a  negative  quantity  for 
ships  with  positive  static  stability. 

The  formulation  is  further  simplified  with  the  introduction  of  the 
two  quantities 


K' 

SB 

*s 

s 

Z' 

w 

X' 

S 

-  \ 

s 

\ 

K'   may  be  regarded  as  a  measure  of  the  control  effectiveness  or  as  the 

ratio  of  the  lift  force  generated  by  the  stern  planes  to  the  lift  force 

generated  by  the  fixed  surfaces  of  the  hull  due  to  angle  of  attack.   X' 

s 

is  the  non-dimensional  effective  radius  at  which  the  normal  force  gen- 
erated by  the  stern  plane  acts.   The  negative  sign  arises  from  the 
coordinate  system  employed. 

Making  all  of  the  above  substitutions  in  equations  (l)  and  (2) 
after  dividing  through  by  m'  yields  the  equations  of  motion  of  the 
hull  with  which  we  shall  be  concerned. 

(3)  (1  +  k2)  w'  -  Zk_w«  -  (  1  +  x'  Z±)   q'  =  K^  Z;   g 

(4)  -d  -  *{   z;  )w'  +  (k^  )2  (1 

m'  m'      u 

The  utility  of  this  formulation  lies  in  the  fact  that  U  is  the 
only  speed-dependent  parameter  of  importance,  and  the  number  of  geome- 
try-dependent parameters  has  been  greatly  reduced.   If  it  is  assumed 

that  Z   varies  as  the  first  power  of  the  speed,  which  is  not  unreason- 
w 


m' 

m'      s 

k')    q' 

"   (X3')2zw     *' 
m« 

+  0 

u 
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able,  then  it  can  be  seen  that  Z  '/m'  is  independent  of  speed  from  the 
following:  w 


Z  ,    ■ =  -  W 

w'     — 

-irr         Urn 

For  hulls  of  normal  design,  the  ranges  of  variation  of  k  ,  k' ,  and 
k'   are  very  limited.   Similarly,  the  effective  moment  arms  and  radii 
X. ' ,  x  ' ,  and  (*,' )   are  confined  to  sufficiently  limited  ranges  of 
values  that  they  may  be  considered  to  be  of  secondary  importance  in 
this  investigation.   Accordingly,  these  parameters  are  assigned  the 
following  normal  values  and  henceforth  treated  as  constants. 

x  •  -  -0.5         (k  ')2  -  °-°45 

x2'  =  0.5  k2    .  0.933 

(x  ')2  =  0.25  k'    -  0.805 

x  •  =  -0.5 
s 

Then  the  dynamic  characteristics  of  the  hulls  to  be  considered 

«^      , 
depend  upon  the  three  primary  parameters  U,   Zw'/m'»   and  K  '.   As 

0 

noted  above,  the  metacentric  height  Z   also  enters  implicitly  as  a 
forth  major  parameter.   These  primary  parameters  will  be  considered  to 
fluctuate  about  assigned  normal  values  as  follows: 


u 

S3 

5.0  i 

4.0 

Zw' 

as 

-1.5  i 

0„5 

m' 

K' 
s 

33 

0.3  ± 

0.1 

It  is  considered  that  the  indicated  ranges  of  Zw'/m1  and  K 


j 


3 


include  most  body  of  revolution  hulls  existing  e)r  to  be  constructed. 
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The  actual  speed  range  indicated  by  the  range  of  U  depends,  of  course, 
on  the  value  of  metacentric  height  selected.   This  figure  normally  lies 
in  the  range  -0.5  to  -2.0  feet.   For  Z_.  =  -0.5  feet,  the  indicated  range 

.D 

of  U  corresponds  to  a  speed  range  of  2.3  to  21  knots.   For  Z_  ■  -2.0 
feet,  this  range  corresponds  to  speeds  from  4.7  to  41  knots.   Again,  it 
is  considered  that  the  ran^e  0f  U  considered  includes  all  combina- 
tions of  Z_,  and  U  of  practical  significance. 

a 

C.  Formulation  of  the  Proposed  System 

By  applying  Laplace  transform  operations,  in  which  the  operator  "s" 
implies  differentiation  with  respect  to  non-dimensional  time   ,  to 
equations  (3)  and  (4)»  one  obtains  the  following  transfer  functions  for 
the  hull. 
(5)  6(s)   -    (E.s  +  F  )  s 


AsJ   +  Bs+Cs  +  D 


(6)        y    (b)        .       Nx   s2  -  K2   s  -  N3 

§s   ^B'  As3   +  Bs2   +   Cs  +   D 

where 

A  =      (l+k2)(ky«)2   (l+k«) 


B  =   -Z 


Wl   f(k'y)2  (l+k«)  +  (l+k2)   (x3')2j 


(x    ,-x1«)    +   1   +  k  -   1 


^2 
m'  IT 


-   Z    '    1 
w 

m*      -2 

-(l+k2)  x( 

1    K    1 

3        s 

Z    ' 
w 

m' 

Z   '      K   ' 
w          s 

"  Z    ' 
w 

m' 

<V 

m' 

F  = 
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n2 


Nl   =  V  V   (V}   (l+k'} 


m' 


N, 


K 


(x  ■  )  +xs'  x2«   -k2xs' 


N,   -  K  •  Z  ' 
3      s   w 


•  Z  •  f  Z  •  r 

w  J   w    I 

m'  1     m'    L 

(x  •  -  x  «)Z  •  -  1   +  ll 

s     1  _w_ 

m«    TT2    J 


Equations  (5)  and  (6)  are  rewritten  in  standard  from  as  follows: 
(7)   g  (s)      K^s+fOs 


S   (s 

vsv 


(8)   Z  •  (s) 

6 


3  t  2  _ 
s  +  J^s  +  J,  s  +  J 
2     1     o 


K2  (a  +  E2s+  P2) 

3  r  2    ;       7 

s  +  J^s  +  J„  s+  J 


where 


K..  ■  E_ 

IT 

N  =  F 

~e7 


K2 

=  N1 

A 

E2 

--"2 

N. 


J2  -_B 
A 

J,   -  C 


2   =  -N. 


N, 


J    -  B 
o 


Substituting  the  assigned  values  of  the  secondary  parameters 
given  above,  these  coefficients  are  expressed  as  functions  of  the 
primary  parameters  as  follows: 


Kn  -  6.16  K  *  Z  ' 
1         s   w 

m' 


K2  =  0.517  Ka'  Zw« 


n;' 
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N  =     1.034 

E2   = 

J0=*  -3.6  Z  ' 

c.                         W 

F2   = 

m' 

Jx-     6.37^  l.?33  - 

Z   •   -1^ 
w        y 

J      = 
o 

U2 

m» 

-5.745 
-12.3  (1-1) 

a2 


-6.37  V 


^2 
U 


m' 


The  following  points  are  now  evident  concerning  the  transfer 
functions  defined  by  equations  (7)  and  (8): 

1.  both  gains  are  always  negative  since  Zw1  is  a  negative 

m' 
number  and  K  '  is  always  positive. 

2.  both  are  stable  transfer  functions  in  that  neither  has 
a  right  half  plane  pole. 

3.  (8)  is  non- minimum  phase  since  E  is  negative. 
Using  these  relations,  the  system  to  be  investigated  is  re- 
presented by  the  block  diagram  of  Pig.  IV. 


*-*Q 


(9) 


Pig.  IV 

The  open-loop  relation  for  the  unity  feedback  system  is: 

K101  (s) 


bi 


(b) 


(b) 


K  .  G 
cl  c 


1  (s)   i-Wi°o2  ra   v8  <°> 
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where   K   .0        and  K   _Q        are    the   compensation  functions    to   be   formu- 
cl   ci  c2   c2 

lated,    and 

K1G]_   (s)      =  K]_   (s+W)    s 


3      t       2      T  _ 

s  +   J_    s  +   J.    S+   J 
2  1  o 

K2G2   (s)      =       K2   (s2  +  E2   s  +  P2) 


EI    (s  +   N)    s 


e(s)      .        Z'  .    (s)   -  Z'      (s) 

gi  go        ' 

Substituting  the  expressions  for  K..    G..    ( s)   and  K     G     (s)   into 
(9),   one^gfbtains 

(10)        ■•        (s)      _  KolK20ol   (s)    (s2+E2s+P2) 


3  ? 

e    (s)  b     +  J.  ■     +  J,    i  +  J     -  L   I  ,  8     0(s)(s+N)s 

2  1  o  1      c2        c2 

Now  we  choose  G  _  (s)  to  be  of  the  form 

C  i— 

G„(s)=    s  +  as  +  b 
02        s  (s  +  N) 

in  which  we  note  that  N  is  a  constant  and  therefore  independent  of 

parameter  variations. 

With  G   (s)  of  this  form,  (10)  becomes 

C  c. 

(11)   i-   (s)   =  K2  I       0ol  (s)  (s2+  E23  ♦  P2) 


^o 


e  W      s3+  (J2-  KlK<j2)s2+  (Jr  KlKo2a)S  +(Jo-KlKo2b) 

and  the  utility  of  this  selection  of  G  0(s)  is  immediately  obvious.   By 
suitably  choosing  K  _,  a,  and  b,  the  poles  of  the  open-loop  relation 

C  c. 

may  be  conveniently  located  in  the  s  plane.   (As  will  be  indicated 
later,  the  situation  is  not  quite  so  straightforward  due  to  saturation 
inSa.) 
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As  stated  above,  the  system  will  "be  designed  with  the  primary 
parameters  assigned  their  normal  or  mid-range  values.   Inserting  these 
numbers  in  (ll)  yields 

(12)   Z'  (s)   -•233KolGcl(e)(s2-  5.745s  -  11.8) 

go    -  

•  e  (a)   s3  +(5.4  +  2.77K  J  s2  +  (3.68  +  2.77K  .a)  s  + 

Cd  Cd 

(.382  +  2.77  Kc2b) 
For  the  present,  take  G  (s)  equal  to  unity.   Then  the  zeros 
of  the  open-loop  relations  are  located  at  +  7.35  and  -  1.61.   Since 
K   must  be  positive,  we  are  concerned  with  variations  of  the  negative 

GJL 

gain  factor  -.233  K  ..  as  far  as  stability  and  performance  of  the  closed- 
loop  system  are  concerned.   The  locations  of  the  poles  of  (12)  are 
determined  by  experimenting  with  the  root  locus  plot  of  this  relation 
to  achieve  the  desired  resonant  frequency  and  damping  ratio  of  the 
dominant  pole-pair  of  the  closed-loop  function.   The  resulting  config- 
uration is  shown  in  Pig.  V.   It  is  to  be  noted  that  the  most  negative 
pole  of  (12)  is  one  of  the  dominant  poles. 

With  the  poles  of  (12)  located  at  -1,  -4,  and  -9  as  shown,  the 
denominator  of  (12)  becomes 

s3+  14s2+  49s  +  36  -  s3+  (5.4  +  2.77Kq2)s2  +  (3.68  +  2.77  Kc2a)s  + 

(.382  +  2o77Kc2b) 

Equating  coefficients  of  the  various  powers  of  s  yields 

K  r.   «■  3.1 

C<£ 

a   -  5.26 
b   =  4.14 

The  phase-magnitude  plot  of  the  open-loop  relation  incorporating 
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these  values  is  shown  in  Fig.  VI  as  curve  A.   From  this  plot  it  is 

evident  that  G  .  should  be  a  lag  network  in  order  to  increase  the  low 

frequency  gain  of  the  system.   Taking  G  ,  (s)  to  be  given  by 

cl 


G  ,  (s)  «  s 

cl       — 


+  2 


s  +  .2 

curve  A  is  modified  as  shown  in  curve  B  of  Fig.  VI.   Applying  curve 

B  to  a  Nichol's  chart  indicates  a  required  gain  factor  of  14  db  for 

an  M  of  2.3  db.   That  is 
P 
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log  .233  Kcl  =  14 


K  •  =  21.4 
c 

The  frequency  response  of  the  closed-loop  system  is  shown  in  Fig.  VII, 

A  steady  state  error  of  -.4  db,  or  4.5$>  is  indicated  in  depth  rate, 

but  this  is  considered  acceptable. 

With  K  ..  G  ,  (s)  and  K  0  G  _  (s)  thus  determined,  the  overall 
cl  ol         c2  c2  * 

system  was  simulated  on  a  Reeves  Electronis  Analog  Computer.   As  will 

be  indicated,  it  was  necessary  to  modify  0  .  (s)  in  order  to  realize 

a  satisfactory  step  response.   The  computer  set-up  is  described  in 

Appendix  A,  and  the  plotted  responses  in  §  ,0,  and  e  for  a  step 

input  of  Z' .  are  shown  in  Appendix  B. 
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Ill 

RESULTS 

A.  General 

The  proposed  system  was  simulated  on  a  Reeves  Electronic  Analog 
Computer  as  described  in  Appendix  A.   A  constraint  of  plus  and  minus 
25  degrees  or  .436  radians  was  placed  on  the  maximum  stern  plane  angle 
permissible.   The  input  <<fco  the  system  was  chosen  to  be  a  step  function. 

Z'±  -  .5  u.^r) 

This  choice  was  made  on  the  basis  that  in  a  steady  glide  w'  should  be 
desirably  zero  making 

• 

z'  =  -e 

q 

where  0  is  the  maximum  permissible  steady  pitch  angle.   This  value  of 
0  was  taken  to  be  -.5  radian  or  -28.6  .   30  is  usually  taken  as  the 
maximum  allowable  value  of  0  in  medium  speed  submarines.   The  actual 
peak  value  of  0  indicated  in  the  results  of  Appendix  B  are  only  ap- 
proximate due  to  the  assumption  that  0  =  sin  0.   The  error  introduced 
by  this  approximation  is  justified  in  tie  sense  that  using  a  high 
value  of  0  permitted  observing  the  system's  general  performance  in  the 
most  extreme  condition  normally  encountered. 

B.  Evaluation  of  Results 

With  the  system  as  designed  in  Chapter  II  set  up  on  the  computer, 
and  the  primary  parameters  assigned  their  normal  values,  it  was  observed 
that  the  overshoot  in  0  and  e  were  of  the  order  of  thirty  percent.   This 
was  due  to  excessive  lag  introduced  by  G   and  saturation  in  £*  .   Ac- 

Ci  S 

cordingly,  the  locations  of  the  pole-zero  pair  of  the  lag  network  were 
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shifted  to  -0.1  and  -1.0  respectively  by  changing  the  setting  of  po- 
tentiometer number  15.   This  corresponded  to  changing  the  value  of 
k  in  the  lag  function  representation  given  in  Appendix  A  from  1  to  .5. 
The  effect  of  this  adjustment  on  the  phase-magnitude  plot  of  the  open- 
loop  function  and  on  the  closed-loop  response  are  shown  in  Curve  C  of 
Pigs.  VI  and  VII.   After  this  adjustment,  the  system  performed  satis- 
factorily.  Plotted  results  are  shown  in  Appendix  Bc 

The  following  observations  are  to  be  noted  concerning  these 
results: 

1.  the  system  specifications  are  met  for  values  of  Z  ,    of 

-1.5  and  -2.0  for  all  values  of  U  greater  than  2.0  and 
all  values  of  K^  considered. 

2.  the  system  does  not  function  satisfactorily  for  Z  '  =-1.0. 

"  w 

m' 
At  values  of  U  of  1.0  and  2.0,  the  system  either  responds  very 
slowly  and  with  large  steady-state  error,  or  shows  very  little  res- 
ponse at  all,  especially  for  low  values  of  K' .   This  is  due  directly 
to  the  small  lift  force  generated  by  the  stern  planes  at  slow  speed 
and  cannot  be  taken  as  a  failure  of  the  plane  angle  control  system. 
At  these  speeds,  S     takes  on  a  maximum  value  and  retains  this  value 
as  would  be  expected.   For  this  reason  values  of  U  less  than  3.0  are 
not  considered  in  the  results. 
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IV 

DISCUSSION  OF  RESULTS 

It  was  noted  that  the  system  did  not  function  satisfactorily 

for  Z'/m'  equal  to  -1.0.   The  reason  for  this  becomes  apparent 

when  the  effects  of  this  parameter  and  of  saturation  is  5     on  the 

locations  of  the  open-loop  poles  of  the  system  are  considered. 

Limiting  the  maximum  value  of  §  is  equivalent  to  limiting  the  ef- 

s 

fectiveness  of  the  compensation  gains  K  n  and  K  „.   It  is  recalled 

°      cl      c2 

that  the  open-loop  poles  of  the  system  were  moved  further  out  into 

the  left  half  of  the  s  p^ane  by  a  suitable  choice  of  K   .   The  effect 

of  saturation  in  5  is, therefore,  to  constrain  these  poles  to  lie  to 

s 

the  right  of  their  desired  locations.   The  effect  is  aggravated  by 

the  tendency  of  J,  as  well  as  J  to  approach  zero  as  U  increases  when 

Z  'An'  equals  -1.0.   For  values  of  Z  ,/ml  between  zero  and  -1.0,  J, 

approaches  a  negative  number.   In  the  extreme  case,  two  of  the  poles 

fall  between  the  origin  and  the  zero  of  the  lag  network,  and  these 

become  the  dominant  pole  pair.   The  net  result-is  a  lowering  of  the 

closed-loop  natural  frequency  and  damping  for  a  given  gain.   Two 

possible  remedies  present  themselves.   First,  the  pre-saturation  gain 

K  ..  may  be  decreased  to  keep  the  system  in  the  linear  region  of  8 

longer,  with  a  resultant  decrease  in  steady-state  accuracy.   Fig.  XXVI 

of  Appendix  B  shows  this  result  for  K  n=  9*25  instead  of  21.4.   A 

second  approach  lies  in  increasing  the  post-saturation  gain  of  the 

hull  via  K' .   This  is  equivalent  to  supplying  more  control  surface 
s 

area.   Figs.  VIII  through  XIII  of  Appendix  B  indicate  this  tendency. 
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A  more  detailed  discussion  of  the  effects  of  saturation  indicated 
would  require  a  non-linear  analysis  of  the  system  and  will  not  be 
attempted  here. 

The  question  immediately  arises  as  to  why  these  effects  were 
not  predominant  at  the  higher  values  of  Z  '  investigated.   The 
necessity  of  shifting  the  zero  of  the  lag  network  to  the  right  men- 
tioned above  indicates  that  these  effects  were  present  when  Z  '/m' 

w  ' 

equalled  -1.5,  but  a  suitable  remedy  was  quite  simply  achieved.   The 

effect  is  not  so  pronounced  simply  because  J.,  and  J  assume  larger 

values  and  the  hull  gain  is  larger  for  these  values  of  Z  ,/m'» 

This  failure  of  the  system  at  low  values  of  Z  '  is  not  considered 

w 

a  serious  limitation  on  the  basic  approach  presented.   For  normal  hull 

/  1*5 

forms,  the  magnitude  of  Z  ,/m9  is  rarely  less  than  1.2.     The  re- 
sults simply  confirm  the  known  fact  that  hulls  with  little  directional 
stability  must  be  provided  with  larger  control  surfaces. 

The  remainder  of  the  plotted  results  of  Appendix  B  for  values 
of  Z  '/m'  of  -1.5  and  -2.0  are  considered  to  justify  the  method  pre- 
sented.  The  behavior  of  pitch  angle,  a  critical  quantity,  in  the  dive 
entry  maneuver  simulated  is  considered  to  be  quite  satisfactory.   The 
error  in  steady  state  depth  rate  is  maintained  less  than  eight  percent 
in  all  instances,  and  the  overshoot  is  negligible.   It  will  be  re- 
called that  the  calculated  steady  state  error  of  the  system  as  com- 
pensated was  four  and  one-half  percent.   From  Fig.  XVII,  the  steady 
state  error  of  the  system  as  designed  is  seen  to  be  four  percent.   Since 

U  appears  in  the  denominator  of  the  expressions  for  J  and  J  ,  decreas- 

V 

ing  U  will  move  the  poles  of  the  open-loop  relations  to  the  lef o  and 
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thus  increase  the  steady  state  error,  an  effect  clearly  shown  in  the 
plotted  results.  The   converse  is  true  when  U  is  increased.   It  is 
seen  from  the  plots  that  this  effect  is  made  quite  small  hy  the 
compensation  technique  employed. 

It  will  be  recalled  that  the  requirement  for  a  steady-state 
plane  angle  to  counteract  the  metacentric  righting  moment  was  neg- 
lected in  the  system  formulation.   From  the  results,  it  is  seen  that 
this  requirement  is  almost  completely  satisfied  by  the  system  due  to 
the  effect  of  metacentric  righting  moment  on  depth  rate.   For  low 
values  of  U,  an  appreciable  steady  state  plane  angle  exists  but 
there  is  still  a  small  pitch  rate. 

It  is  noted  that  the  effect  on  steady  state  error  due  to  changes 
in  hull  gain  through  variations  in  K'  and  Z  '/m'  are  smaller  than  the 
effect  of  metacentric  moment  noted  above.  However,  several  runs  were 
made  with  K  ^    adjusted  to  correct  for  this  change  in  hull  gain.   The 

C-L 

runs  were  made  only  for  the  case  U  ■  5»0>  but  "the  results  shown  in 

Figs.  XXVII  and  XXVIII  indicate  that  a  slight  improvement  in  steady 

state  error  results. 

In  the  basic  system  presented,  it  was  assumed  that  an  infinite 

plane  rate,£  ,  was  available.   This  is  not  the  case  in  practice.   The 

compensation  method  presented,  however,  can  easily  be  extended  to 

handle  the  case  when  plane  actuator  dynamics  are  included  in  the  K.G 

function.   For  example,  if  these  dynamics  are  representable  by  a 

transfer  function  of  the  form  ,  addition  of  the  term  (p  s+l) 

p  s+l 

in  the  numerator  of  both  G  .,  and  G  _  would  cancel  this  function  and 

cl      c2 

require  minor  changes  in  the  constants  associated  with  K  ,G  .,  and 

°  cl  cl 
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K  Jj      .   Again,  complete  cancellation  would  be  achieved  only  in  the 
linear  region  of  stern  plane  operation. 

It  is  becoming  common  practice  to  evaluate  performance  of  a 
proposed  submarine  control  system  in  terms  of  the  comparison  between 
the  trajectories  achieved  for  step  inputs  to  the  system  and  experi- 
mentally determined  minimum  time  or  limit  trajectories.   Due  to  time 
and  computer  capacity  limitations,  such  a  direct  comparison  was  not 
carried  out  in  this  investigation.   It  is  maintained,  however,  that 
on  the  basis  of  the  indicated  behavior  of  stern  plane  angle  the 
plotted  trajectories  of  depth  rate  error  and  pitch  angle  for  the 
cases  of  interest  closely  approach  those  observed  in  the  limit  dive 
entry.   In  the  results  of  Appendix  B  for  the  cases  in  which  Z  '/m1 
assumes  the  values  of  -1.5  and  -2.0,  S      shows  at  most  one  complete 
reversal  from  its  initial  maximum  value  before  it  quickly  settles 
out  to  its  steady  state  value.   This  behavior  is  necessarily  charac- 
teristic of  a  minimum  time  trajectory.   More  will  be  said  on  this 
point  in  the  following  chapter. 
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A  PEOPOSED  PULL-OUT  TIME  PREDICTOR 

The  proposed  system  functions  satisfactorily  in  carrying  the  ship 
through  a  dive  or  rise  entry  phase  and  maintaining  the  required  dive 
or  rise,  but  it  cannot  inform  a  human  operator  when  to  pull  out  to 
achieve  a  desired  depth.   As  was  indicated,  a  computer  or  programmer 
could  be  substituted  for  the  human  operator,  but  a  simpler  solution 
would  be  desirable.   In  the  course  of  the  investigation  reported 
herein,  a  simple  analog  method  of  computing  pull-out  time  presented 
itself.   The  basic  idea  stems  from  an  optimum  controller  concept 

proposed  by  Lieutenant  M.  C.  Davis  USN  in  his  doctoral  thesis  on  multi- 

17 
dimensional  random  processes.    This  technique  is  considered  worthy 

of  note  and  further  investigation  as  an  auxiliary  component  of  the 
proposed  system. 

If  the  optimum  trajectory  in  depth  is  specified  to  be  one  executed 
in  minimum  time  with  no  overshoot,  then  the  optimum  controller  applies 
maximum  control  effort  throughout  the  entry  and  pull-out  phases. 
Consider  a  ship  controlled  by  the  proposed  system  to  be  in  some  steady 
dive.   In  pull-out,  the  operator  could  simply  return  his  control  col- 
umn to  the  zero  position,  commanding  zero  depth  rate,  or  he  could  com- 
mand varying  amounts  of  negative  depth  rate  to  speed  up  the  pull-out. 
But  we  have  seen  that  the  system  presented  in  this  paper  approaches 
the  performance  of  the  optimum  controller  itself  when  the  input  is  the 
desired  step  change  in  depth  rate.   Therefore,  commanding  a  negative 
depth  rate  in  pull-out  cannot  give  a  faster  response  than  simply 
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returning  the  control  column  to  zero.   The  problem,  then,  is  deter- 
mining when  the  operator  should  zero  his  control  column  to  level  out 
at  a  desired  depth. 

Th  problem  and  an  indication  of  a  possible  solution  are  shown  in 
the  following  figures.   In  (a),  the  ship  is  shown  in  some  steady  dive, 


iqo 


zero 
cfetectbr 


O 


//rcS/Cd/vr 


as  indicated  by  the  solid  line,  which  will  carry  it  through  the  de- 
sired depth  Z  ,  .   At  successive  instants,  t. ,  the  present  state  of 
gi  x 

the  system cas  determined  by  speed,  depth,  depth  rate,  pitch,  and 
pitch  rate  are  introduced  as  initial  conditions  into  a  fast  time 
analog  of  the  system,  shown  in  (b).   The  input  to  the  analog  is  a 
step  function  of  equal  magnitude  but  opposite  sign  to  the  existing 
commanded  depth  rate,  a  step  function  returning  commanded  depth 
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rate  to  zero.   The  computing  periods,  in  real  time,  are  short  relative 
to  the  response  time  of  the  hull.   The  output  of  the  analog  is  inte- 
grated to  yield,  in  fast  time,  the  trajectory  in  depth  which  would 
result  if  the  control  column  were  zeroed  at  the  associated  sampling 
instant  t. .   This  projected  depth  response  is  subtracted  from  the 
desired  depth,  represented  by  a  constant  quantity,  and  the  difference 
is  supplied  to  a  zero-crossing  detector.   At  the  first  zero-crossing 
of  this  difference,  the  indicator  is  actuated.   Now  it  is  noted  from 
the  dashed  representations  of  the  projected  trajectories  in  (a)  that 
at  the  first  zero  in  this  difference  both  the  depth  error  and  its 
first  derivative  are  zero.   That  this  will  be  the  case  is  indicated 
by  the  plots  of  depth  rate  error  response  in  Appendix  B.   This  simp- 
lifies the  operator's  task  since  if  he  zeroes  his  control  column  at 
the  first  indication  of  a  zero  crossing  he  will  come  out  on  depth 
with  depth  rate  passing  through  zero.   Then  he  need  only  correct  for 
the  error  in  depth  due  to  the  small  overshoot  of  the  depth  rate  con- 
trol system.   Note  that  pullout  time  is  based  on  the  projected  depth 
at  peak  overshoot  in  depth  rate  so  that  in  a  dive  this  error  will 
always  be  on  the  safe  side. 
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VI 
CONCLUSIONS  AND  RECOMMENDATIONS 

On  the  basis  of  the  results  shown  in  Appendix  B,  it  is  concluded 
that  the  system  presented  is  entirely  feasible  and  worthy  of  further 
investigation  and  implementation.   It  is  noted  that  this  system  pos- 
sesses the  following  very  desirable  characteristics  from  the  operator's 
viewpoint: 

1.  the  system  adjusts  itself  to  any  steady  dive,  or  rise 
commanded  by  the  operator  in  a  minimum  time. 

2.  response  becomes  faster  as  speed  increases,  but  does 
not  become  oscillatory  or  unstable. 

3.  pitch  angle  takes  care  of  itself  within  the  system, 
leaving  the  operator  concerned  only  with  the  primary 
quantity,  depth  rate. 

4.  the  control  system  minimizes  the  effects  of  speed  and 
hull  geometry  variations  by  increasing  the  natural 
frequency  of  the  control  loop  by  a  factor  of  at  least 
five. 

It  is  considered  that  this  system,  in  conjunction  with  a  pull-out 
time  indicator  discussed  in  the  proceeding  chapter,  will  provide  near 
optimum  depth  changing  and  depth-keeping  capability.   Due  to  its 
broadband  nature,  one  system  of  this  type  will  provide  adequate  con- 
trol characteristics  in  a  variety  of  hulls,  obviating  the  need  to 
tailor-design  a  control  system  for  each  new  ship. 
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APPENDIX  A 

COMPUTER  STUDY 
A.   Computer  Simulation  of  the  Proposed  System 

The  system  as  shown  in  Pig.  IV  was  simulated  on  the  analog 
computer  in  the  following  manner.   To  avoid  loosing  sight  of  the 
basic  problem,  the  force  and  moment  equations  of  the  hydroframe  as 
given  in  equations  (3)  and  (4)  were  used  to  represent  the  hull.   To 
facilitate  the  computer  solution,  these  equations  were  rewritten  in 
the  following  form 

(3a)   100  |~  -  .517  j  100  ^°_  -  10   . ll  Zw'   [.5  (100  dO)  +  K^  (10 

L      dr         \-^T      L       "3 — 

m'   u 

>]} 

(4a)   iocwf2_=  -i23(100w,  _  logo  -10  .ll  VI  1 .25  (1000  *-j  ♦ 


+  .i(ioow')1 


dr 


2        ~~~"     TT2  ^     "  -rr    '  v     dr 


dr        L         U 


m' 


d2 

,2 
-      100  £-__ 
dr 

dt2 

+  .5  (ioow')  +  .5  k;  (io)(io  ) 

where     r=  10  *^ 

d  -  10  I — 

Tb  dr 

A  point  concerning  the  analog  treatment  of  the  lag  network  is 
worthy  of  note.  Writing  the  defining  relations  of  this  component  in 
the  following  form, 

gal  (s)  -  Kcl  <B  +  2K) 
e  (s)        s  +  ,2K 

and  carrying  through  the  indicated  operations  and  change  of  time  base 
yields 

10  dgsl   =  K  .  (10  —\   +  K  (2K  ne  -  .2   .) 

— -, —     cl      dr  '        cl       si' 
dr 
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de 
Integrating  once  since  —  is  unavailable,  and  noting  that  the 

initial  conditions  on  the  integration  are  zero  yields 

10  S  .  -  K  ,   (10  e)  +  K  )r2K  .  e  -  .2  §  .1  dr 
si    ol  J  L  ol        ^slj 

or  £sl  -  .olKcl   (100  e)  +  .IK  r[".o2Kcl  (100  e)  -  .2<Sldr 

Now  the  coefficient  of  the  integral  will  appear  in  the  computer 
diagram  simply  as  a  potentiometer  for  values  of  1(  between  0  and  10. 
K  determines  the  location  of  the  pole-zero  pair  of  the  lag  network 
along  the  negative  real  axis  of  the  s  plane  relative  to  the  selected 
position  at  -2  and  -.2,  the  ratio  of  their  distances  from  the  origin 
remaining  constant.   This  facilitates  a  simple  method  of  experimen- 
tally varying  the  phase  shift  introduced  by  the  lag  network  without 
affecting  the  steady  state  error  of  the  system. 

The  feedback  compensation  is  represented  by  the  following: 

10  d5s2  -  .310  (1000  d26  )  +  .163  (1000  d6_  )  +  .128  ( 100  0)  - 

dr  ,2  dr 

dr 

-.104  (I0§s2  ) 

A  limiter  is  employed  to  restrict  the  maximum  value  of  §  =&-,  +  £  0 

to  +  25  or  .436  radian. 

Z'    is  computed  from  the  relation 
go 

o 

100  z»   -  100  w'  -  100  e 

go 

The  complete  computer  diagram  is  shown  on  Plate  I. 
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APPENDIX  B 
ANALOG  COMPUTER  RESULTS 

Pigs.  VIII  through  XXV  show  the  results  of  the  analog  computer 

study  of  the  proposed  system  as  plots  of  &   ,   8,   and  e  vs  the  non- 

s 

dimensional  time  parameter  2^  for  various  values  of  the  three  primary 
parameters  K' ,  ZV/m' ,  and  U.   The  values  of  these  parameters  con- 
sidered are  as  follows: 

U     3.0,  5.0,  7.0,  9.0 


K  •   0.2,  0.3,  0.4 


V   -1.0,  -1.5,  -2.0 


Curves  for  each  value  of  U  are  shown  on  one  plot  for  each  pair  of 

the  parameters  K'   and  Z  '    considered. 

s        w 

a' 

Fig.  XXVI  shows  the  performance  of  the  system  in  pitch  angle 

and  depth  rate  for  Z  '/mt  "  -1.0,  U  -  5.0,  K  ■  -  0.3,  and  K   =  9.25 
Pigs.  XXVII  and  XXVIII  show  this  data  for  Z  '  of  -1.5  and  -2.0 

D  W 

m1 

and  all  values  of  K  '  considered  with  K  n   adjusted  to  maintain  the 

s  ol    " 

product  K  n  K0  constant* 
el  d. 


32 


\ 


1 


L 


r 


+ 


UO 


<r> 


I 


i.  .  . . >  I .  .  .  ] 

*  ■■'■      ■■■'-    ; 


Q) 


M 
H 


CD 


VO' 


1 

i 

1 

i  s 

1        II 

4j  "?    o  < 

it       H        «■      i 

^ 


N 


QJ 


o  ^ 


X 

5 


03 


. 


OJ 

o 

1 

X 

I 


6 


o 
*6 


#1    ^     *    # 


•«Q 


-K 


s 


VO 


■^ 


i 


c\* 


V>      oj      o    0    o 

f?     d      oj    k>  (S 


ii       H 


n     ii 


fx 


||  V    ^  O  O 


-co 


\ 


■-tf 


CD 


X 


1 

/ 
/ 
/ 

\ 


0     0    0 

o6  *o  cs 


{?     1       o     o     o 


II 


It 


^U     V        ^    ^ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


A 


/ 


/A 

//y 

■**- 


|-<r> 


-<vj 


^ 


-N 


S 
X 

1 


<o 


<^ 


o 

CD 


1 

i 

1 
0      ( 

■ 

> 

3    0 

1 

ii 

li        ii 

i 

£ 

* 


V 


uOv 


:vnt\htbrt$T Frtftf f;H Hrm ?  - 


...    •;-■-    [ 

\          > 

o     o    o 

^       0 

\ 

ll 

0)      V>*    <S 

•i        ii      ii                  i 

m  g  |  ?^  |       i 

V 


^ 


£ 

V 

0) 


o 


i 

1 

9    ^      o    c 

)     0 

i     " 

05     6        oi    V, 

' 

i 

II           "              H          1 

m  ^  s ,: 

' 

ijffitftlS- 

- 

/ 

*o 


/ 


/ 


/ 


/ 


^ 


*i 


CD 

ttixiit.il' 


00 


^ 


v© 


V> 


K 


X 

1 


?Z^  c>o 


"0 


iL.'.J  .    _.;..'  .ililiii'L  '— "."i..--ll  =  ::.I-:Iil 


*0 


Q) 


01 

^ 


fO 


?* 


a 


ai 

. 

s 

-w — ^~ -J ;„: :,:, tirrti  v~  .-_: : L. ..^.....,._-..„._^__i_.:.-l_:..;J.  ..-.'.L...uJ._l'..-u-  ...i». 


i 


1 

1 
1 

<vi     o 
i 

m  « 

i 
i 

1            j 
1            1 

o    o    o 
n     ii    " 

^ 


o  ^ 


ft 


\ 


<\ 


o 


I 


»* 


-ro 


** 


•-CU 


X 
X 

i 


\ 

\ 


\ 


CO" 


<D 


JUTITiTTT- 


. 


I 

1 

1 

i 
1 

o 

II 

^     o 

ii 

0 

Ii 

ii      J' 

b 

?^> 

o 


^ 


/ 

•CO 


^o 


<* 


—  o 


>< 

CD 

5! 


CD 


2 


r 


s~ 


2 

4*  r    ai 

0 
3 

u    r    *^< 

0 

IBS 

fep  * 

T 


FIG.    XYV/ 


R* 

to 

03 


cxa 

> 

X 


k 

CO 


I 


a 


ci 

ii 


<v 


l3l?      *" 


o 


o 


■oo 


t\ 


I 


/ 


/ 


s 


A 

/ 

s 

ifi 


^ 


CD 


•> 


' 

* 

si 

_.„. s.^1^ i J L 


APPENDIX  C 
BIBLIOGRAPHY 

1.  "Nomenclature  for  Treating  the  Motion  of  a  Submerged  Body  Through 
a  Fluid,"  The  Society  of  Naval  Architects  and  Marine  Engineers, 
Teohnical  and  Research  Bulletin  No.  1-5  (Apr.  1950) • 

2.  Hoyt,  E.  D.  and  Imlay,  P.  H. ,  "The  Influence  of  Metacentric 
Stability  on  the  Dynamic  Longitudinal  Stability  of  a  Submarine", 
Taylor  Model  Basin  Report  C  -  158  (Unclassified) 

3.  Segalin,  7.  0.,  "An  Analytical  Formulation  of  the  Synthesis 
Problem  of  Corrective  Devices  in  Linear  Servosystems, "  Auto- 
mation and  Remote  Control  (USSR),  vol.  XIX,  No.  2,  Feb.  1958. 

4.  J.  G.  Truxal,  "Automatic  Feedback  Control  System  Synthesis", 
McGraw-Hill  Book  Company,  Inc,  1955» 

5.  Thornbury,  J.  W.  and  Kauffman,  H.  A. ,  "Automatic  Depth  Control 
for  Submarines",  MIT:  M.  S.  Thesis  1949  (Declassified) 

6.  Williamson,  R.  R.  and  Peters,  B.  H. ,  "Trajectory  Analysis  of  the 
Relationship  of  Stability  to  Maneuverability  in  the  Manual  Control 
of  Submarines",  Report  No.  389,  Experimental  Towing  Tank,  Stevens 
Institute  of  Technology,  1949.  (Declassified) 

7.  "Final  Report  Ship  Control  Information  Study  (SUBIC)",  Report  No. 
IM  -  1306  -  V  -  2,  Cornell  Aeronautical  Laboratory  Inc.  (1959). 

8.  "SUBIC  Progress  Broshure,  1  October  I960",  Electric  Boat  Co. 

9.  Pode,  L.  and  Church,  J.  W. ,  "A  Representation  of  the  Submerged 
Dynamic  Behavior  of  Submarines  Proposed  for  the  Special  Devices 
Center  Submarine  Simulator",  DTMB  Report  C  -661  ( CONFIDENTIAL) 


33 


10.  Pode,  L.  and  Church,  J.  W. ,  "A  Study  of  the  Effects  of  Stability 
and  Control  Parameters  on  the  Performance  and  Stability  of  Sub- 
marines in  Pitching  Motion",  DTMB  Report  C  -626  (CONFIDENTIAL) 

11.  Centner,  R.  M.  ,  "Manual  Depth  Control  of  the  AGSS  569  Sub- 
marine", MIT;  M.  S.  Thesis  1957  (CONFIDENTIAL) 

12.  Landweber,  L.  and  Johnson,  J.  L. ,  "Prediction  of  Dynamic 
Stability  Derivatives  of  an  Elongated  Body  of  Revolution", 
DTMB  Report  C  -  359,  19 51.  (CONFIDENTIAL) 

13.  Westbrook,  C.  B.  and  McRues,  D.  T. ,  "Handling  Qualities  and 
Pilot  Dynamics",  Aero/Space  Engineering  vol.  18,  No.  5,  May  1959. 

14.  Chestnut,  H.  and  Mayer,  R.  W. ,  "Servo  Mechanisms  and  Regulating 
System  Design",  Vol.  I,  John  Wiley  and  Sons,  New  York  (1957). 

15.  "Stability  and  Control  Data  for  Current  Submarine  Design",  David 
Taylor  Model  Basin  Miscellaneous  Report  105  (1949)  (CONFIDENTIAL) 

16.  Preston,  J.  S.  "Non-linear  Control  of  a  Saturating  Third  Order 
Servome onanism",  MIT;  M.  S.  Thesis  1954. 

17.  Davis,  M.  C. ,  "Models,  States,  and  Dimensions  in  Random  Processes", 
MIT;  Sc.  D.  Thesis  I96I. 


34 


